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Section  1 

PLASMA  THLORY 


A.  DRIFT  CYCLOTRON  INSTABILITY 

J.  K.  Lee 

No  additions  to  theory  this  quarter. 

B.  PLASMAS  WITH  FIELD  REVERSAL 

Douglas  Harned  (Pto^.  C.  K.  B-ixdictC) 


T . 


A simplified  model  is  being  examined  to  analyze  the  stability  of 
the  m=0  mode  in  long-layer  field  reversed  plasmas.  The  instability  is  char- 
acterized by  a high  energy  ion  beam  propagating  through  a relatively  cold 
background  plasma.  It  is  of  interest  to  see  if  this  instability,  which  is 
observed  in  simulation  plasmas,  is  real  or  numerical  in  nature. 

Using  the  infinite  radius  limit  with  a long  layer,  the  layer 
may  be  unwrapped  into  a slab  (see  Fig.  1).  In  the  slab  model,  the  m=0  mode 
corresponds  to  rigid  motion  of  the  slab  beam  toward  one  wall.  Our  model  has 
been  considered  by  Friedman,*  which  consists  of  a six-cell  beam-plasma  system. 
The  beam  propagates  uniformly  through  a region  one  cell  in  width,  straddling 
the  centerline  of  the  plasma.  (In  the  six-cell  model  this  corresponds  to  the 
line  separating  the  third  and  fourth  cells.)  The  walls  are  treated  as  per- 
fect conductors.  The  only  zero-order  magnetic  field  is  assumed  to  be  the  self- 
field (B^)  of  the  zero-order  beam  current  (J°) . Quasineutrality,  ne  = nj+nb 


A.  Friedman,  Cornell  University,  private  communication. 
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(ng,  n.,  and  are  the  electron,  ion  and  beam  densities,  respectively),  is 
assumed.  The  displacement  current  is  neglected  as  we  are  looking  only  at  low 
frequency  phenomena.  The  electrons  are  taken  as  a perfectly  conducting  fluid. 
The  f in i te-d i fferenced,  linearized  equations  of  motion  for  both  the  beam  and 
the  plasma  are  used  to  obtain  a quadratic  equation  for  the  frequency  u>+iy  . 
There  are  two  possible  modes,  symmetric  and  ant i -symmetr i c with  respect  to 
the  perturbed  electric  field.  When  the  roots  are  found,  using  the  root-solver 
ZANALYT,  the  symmetric  mode  appears  stable  while  the  antisymmetric  mode  has 
positive  growth  rates  for  certain  values  of  the  beam  velocity  and  density.  The 
most  unstable  regime  is  for  beams  having  densities  larger  than  the  plasma  den- 
sity and  velocities  large  relative  to  the  Alfven  velocity.  Another  weaker  area 
of  instability  was  found  for  low  density,  low  velocity  beams. 

In  order  to  prevent  problems  that  might  be  a result  of  the  small 
number  of  cells  used  by  Friedman  in  the  six-cell  case,  we  have  expanded  the 
model  by  constructing  an  algorithm  for  an  arbitrary  number  of  cells.  However, 
the  width  of  the  beam  was  retained  at  one  cell.  A polynomial  of  order  N = 1/2-1 
was  obtained,  where  I is  the  number  of  cells.  This  polynomial  was  then 
solved  numerically,  again  using  ZANALYT.  The  results  for  cases  from  8 to  20 
cells  confirmed  the  results  of  the  six-cell  model.  The  growth  rates  were 
comparable  (on  the  order  of  the  ion-cyclotron  frequency)  and  instability  was 
observed  in  the  same  regions  of  parameter  space  (the  parameters  being  beam 
velocity  and  density). 

The  multicell  model  still  contains  the  weakness  of  having  only  a 
one-cell  wide  beam.  This  problem  has  caused  us  to  seek  an  analytic  solution; 
this  work  is  not  yet  complete.  The  model  is  essentially  the  same  as  that 
described  above.  The  linearized  equations  of  motion  are  used  for  the  beam 
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An  important  point  here  is  that  in  the  models  using  a one-cell  beam  the 
v^(B°/c)  and  terms  do  not  appear  because  the  antisymmetry  of  B° 

made  these  terms  vahish  upon  cell  averaging.  The  electrons  are  still  consi 
dered  as  a perfectly  conducting  fluid. 


R 

E 1 + v 1 x — = 0 ; 
e c 


quasineutrality  is  again  assumed.  When  these  relations  are  combined  with  Max- 
well's equations,  second  order  differential  equations  are  obtained  for  each 
region  in  the  model  for  the  perturbed  y directed  electric  field  . These 
equations  can,  in  principal,  then  be  solved  by  using  boundary  conditions  at 
the  perturbed  boundary.  This  solution  is  not  yet  complete. 


Section  II 
SIMULATION 

A.  QUASI  LINEAR  SIMULATION  OF  A MIRROR  MACHINE 

Dr.  Y.  Matsuda  with  Dr.  H.  L.  Berk  (LLL) 

This  project  is  in  limbo;  there  will  be  no  further  report. 

B.  ONE  DIMENSIONAL  ESI  CODE  FOR  SHEATHS 

Y-J  Chen  (P'lorf.  C.  K.  Birdiatt) 

No  special  work  to  report  this  quarter. 

C.  PARTICLE  TRAJECTORIES  IN  A CUSP  FIELD 

Y-J  Chen  (PAO;$.  C.  K.  Bindbatt) 

No  special  work  to  report  this  quarter. 

D.  SIMULATIONS  OF  DRIFT-CYCLOTRON  INSTABILITY 

J.  K.  Lee  (Pao^.  C.  K.  BindAoM) 

We  have  made  several  si mul at  ions  of  the  drift-cyclotron  instab- 
ility using  EZOHAR.  These  results  are- prel iminary;  only  a general  descrip- 
tion will  be  presented  in  the  following. 

Our  simulations  include  exact  particle  dynamics  of  both  electrons 
and  ions  (i.e.,  the  guiding  center  approximation  of  the  electrons  was  not  used) 
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with  flute-like  (k| * 0)  two-dimensional  electrostatic  perturbations.  Typical 
parameters  are 

number  of  spatial  grids  = 32  * 32 

number  of  particles  (electrons  or  ions)  = 16x16x16  to  32x32x32 

we/wci  = “ce^ci=  VWci)2=mi/me=25’  11,  9 (three  cases  50  far) 

id  At  ^ 0 . 5 
pe 

p ./L  =4^1 
i n 

where  the  above  symbols  are  defined  in  the  last  QPR  (p.  2). 

Simulations  show  several  modes  growing  exponentially  in  time 

up  to  u . t 'v  1 0 . The  comparison  of  this  linear  behavior  with  the  existing 

linear  theories  is  not  complete  at  this  moment,  since  most  of  linear  theories 

use  the  local  approximation  where  p./L«  1 . The  only  nonlocal  theory 

1 1 " 

(known  to  us)  is  that  of  ROOTS,  which  has  a sinusoidal  guiding  center  density 
profile  like  our  mode),  but  which  is  periodic  in  Vn  direction,  where  our 
model  is  unbounded  in  that  direction.  The  linear  growth  rates  and  frequencies 
in  the  simulations  are  found  to  agree  better  with  the  nonlocal  theory  predic- 
tions rather  than  with  those  of  the  local  theory. 

At  saturation,  the  total  electrostatic  field  energy  reaches 
about  2%  of  the  initial  ion  kinetic  energy  for  the  above  mass  ratios.  These 
simulation  saturation  levels  are  now  being  compared  with  a small  number  of 
existing  nonlinear  theories,  and  the  preliminary  results  show  quite  good  com- 
parison. At  and  after  saturation,  the  phase  space  plots  (v  - y)  of  electrons 
and  ions  form  vortex-like  structures  in  accordance  with  the  most  unstable  mode 
and  the  density  profiles  of  both  species  start  to  develop  a stretching  pattern 

(both  of  these  are  similarly  observed  in  the  lower-hybrid-drift  \ .ability 
2 

simulation  ).  The  density  gradient  stretching  was  especially  clearly  seen 
in  movies;  these  have  the  potential  contour  plot  at  the  top  of  the  frame  and 
the  density  profiles  of  electrons  and  ions  superimposed  in  different  colors 
at  the  bottom.  The  movies  also  show  clearly  the  mode  localization  in  x (or  y) 
and  the  mode  propagation  in  the  x-y  plane  (mostly  in  the  ion  diamagnetic  drift 
wave  direction). 

'm.  J.  Gerver,  "ROOTS,  A Dispersion  Equation  Solver",  UCB/ERL  Memorandum  No. 

M77/27  (Oct.  31 1 1976). 

2D.  Wlnske  and  P.  C.  Liewer,  "Particle  Simulation  Studies  of  the  Lower  Hybrid 
Orift  Instability",  Phys.  Fluids  21_,  pp.  1017-1075  (1978). 
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Suction  III 

CODE  DEVELOPMENT  AND  MAINTENANCE 

A.  ESI  CODE:  MODIFICATION  OF  ESI  TO  USE  POST-PROCESSOR  ZED 
Dr.  W.  M.  Fawley 


ZED  is  a post-processor  routine  that  allows  the  user  to  manipu- 
late and  plot  time  histories  of  variables  and  mode  amplitudes  contained  in 
the  output  produced  by  the  plasma  simulation  code  EZOHAR.  The  manipulations, 
such  as  Fourier  transforming  in  time  and  cross-correlation,  can  be  extremely 
useful  in  determining  certain  properties,  e.g.  mode  exponential  growth  rates. 
ZED  has  proved  so  beneficial  in  use  with  EZOHAR  that  we  decided  to  modify 
the  simpler  codes  ESI  and  EMI  to  produce  ZED-readable  output. 

ZED  expects  both  a "STATE"  file  and  "HISTORY"  file  to  exist  on 
disk  when  it  is  run.  The  STATE  file  tells  ZED  the  last  time  step  (IT)  com- 
pleted (actually  the  last  time  step  divisible  by  50),  the  value  (DT)  of  the 
time  step,  and  how  often  (in  time  steps)  each  variable  is  stored  (ISMODE, 
normally  1;  i.e.,  each  va- table  is  stored  each  tine  step).  ZED  will  normally 
ask  the  user  for  his/her  box  number  and  file  identification. 

The  HISTORY  file  contains  the  actual  time  histories.  The  file  is 
split  into  two  parts.  The  first  part  contains  the  time  histories  of  vari- 
ables such  as  the  electrostatic,  kinetic,  and  total  energy  while  the  second 
part  contains  time  histories  of  mode  amplitudes  and  variable  values  at  parti- 
cular grid  points.  The  user  has  no  choice  as  to  what  variables  are  stored  in 
the  first  part,  but  does  instruct  the  code  via  input  parameters  as  to  which 
modes  are  to  be  stored  in  the  second  part.  The  HISTORY  file  replaces  the 
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functions  performed  by  the  subroutine  HISTRY  in  ESI  and  EMI.  However, 
"snapshots"  of  the  potential  and  other  variables  are  still  written  in  DD80 
files.  Thus,  the  user  will  use  TEKPLOT  to  see  the  snapshots  produced  at  a 
particular  tine  during  the  simulation  run. 

Most  changes  made  to  ESI  and  EMl  were  required  to  write  the  STATE 
and  HISTORY  files.  The  new,  modified  codes  are  called  "ESZ"  and  "EMZ".  Their 
respective  source  listings,  "SESZ"  and  "SEMZ"  are  global  files  contained  in 
user  number  1221,  directory  .ZOH  . Both  codes  must  be  compiled  with  CHATR, 
and  the  binary  routine  "BEM1"  (contained  in  EM1LIB)  must  be  present.  The 
user  must  also  specify  TV80LIB,  ORDERLIB,  and  CF76LIB  (this  order  is  crucial 
for  EMZ). 

In  general,  it  has  proved  very  efficient  (in  CPU  time)  to  use 
binary  libraries  when  debugging  and  modifying  (by  a factor  of  3 " 5)  these 
codes.  Please  see  the  July  1,  1977  QPR.  PP-  55_57»  for  information  on  main- 
taining binary  libraries  and  loading  procedures  (using  LOD) . For  convenience 
the  binary  libraries  BESZLIB  and  BEMZLIB  are  also  present  in  directory  .ZOH  , 
user  number  1221.  Such  library  contains  BEM1  in  addition  to  the  ESZ  and 
EMZ  binary  routines.  The  source  libraries  SESZLIB  and  SEMZLIB  also  exist  in 
directory  .ZOH  . It  is  strongly  suggested  that  each  time  a user  creates  a new 
run  file  ESZ  or  EMZ,  a symbol  table  also  be  created.  This  allows  one  to 
debug  the  programs  dynamically  with  the  system  routine  DBCTRL. 

For  users  who  are  interested  in  modifying  their  own  versions  of 
ESI  and  EMl  to  produce  ZED-compatible  output,  SESZ  and  SEMZ  are  written  in  a 
way  to  point  out  the  necessary  deletions  and  additions.  First  of  all,  the 
new  subroutines  MDSET,  TSET,  and  WRTSTATE  are  needed  to  control  10  to  the 
STATE  and  HISTORY  files.  TSET  establishes  the  HISTORY  file  and  writes  (through 
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"TSAVE"  and  RSAVE)  values  of  variables  such  as  total  energy  into  the  first 
part  of  the  HISTORY  file.  MDSET  reads  inputs  as  to  which  modes  are  stored 
(more  on  this  later)  and  writes  (via  MWSAVE  and  MRSAVE)  the  second  part  of 
the  time  history  file.  WRTSTATE  creates  a new  STATE  file  (and  destroys  the 
old  STATE  file)  each  time  it  is  called.  At  the  end  of  a run,  it  closes  the 
HISTORY  file.  ESZ  and  EMZ  have  different  versions  for  each  of  the  programs. 

The  user  also  should  eliminate  the  routine  HISTORY  and  PLTHST 
as  these  are  no  longer  needed.  Elsewhere  through  the  routines  common  to 
ESI  and  ESZ  and  EMI  and  EMZ,  deletions  and  additions  are  flagged  with  comment 
statements  in  the  following  ways: 

All  lines  to  be  deleted  have  "CK“  in  columns  I and  2.  Each  line 
(or  group  of  lines)  is  preceded  by  a line  with  "CADD"  in  columns  1 -A; 
the  line  will  tell  you  the  number  of  lines  to  be  added  to  your  code.  In  gen- 
eral, this  method  has  worked  well  (one  success  out  of  one  try).  If  your 
version  of  ESI  uses  the  PUTT  compiler,  you  must  change  all  format  statements 
with  delimited  output  statements,  to  be  quote  (")  delimited  statements, 

in  order  to  be  compatible  with  the  CHATR  compiler. 

INPUT  changes:  The  routine  MDSET  requests  via  the  INPUT  file 
information  on  which  modes  to  store.  This  information  should  go  after  the 
"namelist"  input,  separated  by  a "STOP"  statement  (this  STOP  replaces  the 
SEND  in  PUTT-  compiled  ESI  INPUT  files).  The  modes  histories  requests  should 
be  in  right-adjusted  (05,  15)  format.  A blank  line  ceases  the  requests.  The 
05  input  is  bbbTV  where  T specifies  the  type  of  information,  i.e.,  value  at 
a particular  grid  point,  and  V specifies  the  actual  variable  stored  (i.e.. 


electric  potential). 
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T = 0 
1 

2 

V = 1 

2 

3 

k 

5 


EMZ 

only 

EMZ 

only 


B 

a 


mode  energy 

complex  mode  amplitude  (2  words  stored) 
value  at  a particular  grid  point 

electric  potential  (or  longitudinal  electrostatic  energy 
if  T = 0) 
not  used 

current  density  in  y-direction  (EMZ  only) 
total  charge  density 
electric  field  in  x direction 
electric  field  in  y direction  T = 2 
magnetic  field  in  z direction  T = 2 
wave  amplitude  F polarization  T=1 
wave  amplitude  F polarization  T=1 


The  15  input  specifies  each  the  grid  point  position  (T  = 2)  or  the  mode  number 
(T  = 0 or  T = 1). 


After  the  blank  line  which  closes  the  mode  history  input  goes 
the  species  1 and  2 input  for  the  routine  INIT.  No  "STOP"  is  needed  after 
the  mode  history  blank  line.  Sample  input  files  for  ESZ  and  EMZ  are  in 
SESZLIB  and  SEMZLIB  under  the  name  INPUT.  The  main  source  routines  in  the 
two  libraries  are  called  MAINESZ  and  MAINEMZ,  respectively. 


Use  wi th  ZED. 


A user  may  use  ZED  on  the  HISTORY  and  STATE  files  as  soon  as  the 
latter  is  first  written  (after  the  50—  time  step).  If  one  is  familiar 
with  use  of  ZED  with  ZOHAR  output,  there  should  be  no  difficulties.  "ZED 
REPORT",  which  summarizes  the  available  commands,  may  be  obtained  from 
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F I LEM  RDS  1234  ZEDREP(LF)END  /TV 

TRIX  AC  / T V 

PR  I NT (<N I P ZEDREP  BOX  NNN  ZEDREPORT>) 

END 

Until  the  "MODES"  command  is  given,  the  user  may  access  time  his- 
tories of  momenta  ( PX 1,2;  PY 1,2),  kinetic  energies  ( KE 1,2),  electrostatic 

energy  (ESE) , and  total  energy  (TE)  in  ESZ  output.  For  EMZ  output,  the 
equivalent  variables  are  PX 1,2;  PY1  ,2;  KE1,2;  TE;  EXE  (longitudinal  electro- 
static energy);  EYE  (electric  field  energy  in  y direction);  BZE  (magnetic 
field  energy  in  z direction;  EYLE,  EYRE  (wave  energy  in  left  and  right  polar- 
izations); and  TEMPI,  I0NTMP  (average  electron  and  ion  temperatures  in  elec- 
tron volts) . 

After  the  user  specifies  "MODES"  to  ZED,  times  histories  of  the 
mode  energies  and  amplitudes  may  be  accessed.  To  do  this,  the  user  specifies 
the  same  code  for  a given  mode  (e.g.  READ  01  01  00),  EXCEPT  he' she  MUST 
include  "00"  after  the  mode  or  grid  point  number.  ZED  expects  a 2Jd  output 
and  thus  KY  or  NY  (which  is  00  in  lid  output)  must  be  specified.  In  addition, 
to  read  a complex  amplitude  (e.g.  PHI(KX)),  an  ampersand  must  be  included 
after  the  "00".  This  instructs  ZED  to  read  the  second  (imaginary)  word  of 
the  complex  amplitude. 

RELIABILITY:  ESZ  has  been  tested  on  a several  sample  problems 
and  appears  to  work  well.  No  physics  was  changed  in  this  code  relative  to 
ESI,  so  the  two  codes  should  give  identical  answers  on  a particular  problem. 
EMZ,  however,  has  been  not  tested  at  all  and  a CAVEAT  EMPT0R  must  be  ex- 
tended. Inaddition,  a (hopefully)  faster  and  better  routine  was  put  in  the 
routine  "ACCEL"  to  compute  electron  temperatures.  The  changes  may  have  bugs. 
Thus,  any  user  who  is  so  interested  is  encouraged  to  test  EMZ  extensively. 

Good  luck. 


I 
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B.  EMI  CODE 

Nothing  special  to  report  this  quarter. 

C.  EZOHAR  DEVELOPMENT 

J.  K.  Lee 

r 

Code  is  running. 


D.  ESI  + EFL  CODE 

J.  K.  Lee  [Plo^.  C.  K.  &Ud&alt) 

A report  on  this  code  may  be  prepared  for  the  Journal  of  Computa- 
tional Physics.  We  may  add  work  on  the  Lagrangian  fluid  version. 


E.  HD  PARTICLE-FLUID  HYBRID  FOR  LOWER  HYBRID  DRIFT  WAVE  INSTABILITY 
Y-J  Chen  ( P/to . C.  K.  Bild&alt) 

W'v 

A Hd  hybrid  model  for  the  lower-hybr  i d-d  r'i  f t instability  has 

b-en  studied.  Using  the  "ghost"  ion  method  reported  by  B.  I.  Cohen  and 

G.  R.  Smith,'  ions  in  the  simulation  are  distributed  in  v , v and  y 

x y ' 

explicitly.  Each  ion  represents  many  ghost  ions  moving  along  identical  tra- 


jectories except  for  t ime- i ndependent  displacements  in  x . The  initial  ghost 

ion  position  xq  is  calculated  as  it  is  at  the  x*0  plane  at  the  present 

Bruce  I.  Cohen  and  Gary  R.  Smith,  "Efficient  Method  for  Local  Simulation  of 
Drift  Wave  Instabilities",  Proc.  8th  Conf.  on  Numerical  Simulation  of  Plasmas, 
Monterey  CA,  June  28-30,  1978. 


MMMMCMMM 
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time  t . Then  the  charge  q and  mass  m of  this  ghost  ion  are  modified 
by  a weighting  function  W(xq)  = n(xQ)/no  . If  the  wave  perturbation  is  lo- 
calized at  the  x = 0 plane,  i.e.. 


,2  .2  -2 

k » k » L , 

y x 


where  L is  the  density  scale  length,  the  local  approximation  in  the  den- 
sity gradient  direction  is  used  and  W(xq)  can  be  expressed  as 


W(xq)  = 


1 - — , x < L 

L ° 


. *q>L  • 


Since  the  instability  is  characterized  by 


CO  . 0)  « U)  , 

c i ce 


-2,2  -2 

a . « k « a , 

i i e ’ 


2 ,2 
a « L , 
e 


the  ions  are  treated  as  unmagnetized  particles.  The  x displacement  of  ions 


x(t)  « x + v (t-O)t  - E t2  . 

ox  Z nij  Xq 


If  the  ion  is  at  the  x = 0 plane  at  time  t , we  get 


These  three  equations  can  be  solved  by  a predictor-corrector  scheme  which 

takes  care  of  the  convective  terms  v (a/3 y ) in  them.  For  example,  Eq.  (8) 

o 

becomes 


n,p  n- 1 n-i  n-1 


At 


2 Ay 


n-i  n-i  n-i  n-i 

. v +v  v -v 

dp  x.^,  x.  , y...  y.  . 

- _£  L±i L±-  - o . V*  LL-.i. 

dx  2 o Ay 


(U) 


and 


n n- 1 

P.-  P. 

_J J 

At 


+ v 


n-i^  n-i 

n,p  n,p  , v +v 
P. ’ , - P.  dp  x. x.  , 
J+l  J 2.  _ J+i ni 


2 Ay 


dx 


n-i  n-i 
v - v 

h+L  hdL 

Ay 


(12) 


where  p.  in  Eq.  (12)  is  the  final  solution.  Unfortunately,  this  method 
might  take  a relatively  long  time  to  run  on  the  computer. 


The  second  way  is  to  bring  in  the  electron  effects  through 
Poisson's  equation,  with  the  electron  susceptibility  which  can  include 

all  the  nonuniformity  information.  This  will  be  presented  later. 
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RJET  (MINI-USERS'  CENTER)  DEVELOPMENT 


S.  Au-Yeung  [PX-0&.  C.  K.  &uid&aJLl\ 


Our  RJET  to  the  MFENCC  (Livermore)  consists  of  (see  Fig.  1)  one 
PDP  11/04  computer  that  contains 


1 KD11-D 
1 KYI  1 -LA 
1 M9301-YB 
1 MS11-JP 
1 DD11-DK 
1 BA11-KK 
I M7850 
1 DL11-WA 


Processor  module  (CPU) 

Operator's  console 
Bootstrap/Terminator  module 
MOS  memory 
Expander  backplane 
Mounting  box 
Parity  control ler 

Serial  line  interface  and  Real-time 
clock 


2 DMC11-AR 
2 DMC1 1-OA 
1-2  DZ-ll-A 

1 BM873-YC 


DOCMP  Microprocessor  module 
Line  unit  module 

Asynchronous  multiplexer  (8  asynchronous 
serial  lines/unit) 

Restant  loader  containing  ROM  with 
bootstrap  program 


1 Dataphone  4800  baud 
1 Channel  Interface 

upto  16  EIA  standard  RS232C  interface  terminals 
1 Terminal  used  as  the  console. 


The  log  in  procedure  is: 

(1)  Turn  on  the  terminal 

(2)  Type  in 

a(or  c for  CRAY^.user  #_suff ix^acct  ^^(ctrl-z)  password 
Although  the  default  baud  rate  is  300,  it  can  set  up  to  9600  by: 


* 


> 

: 


(1)  Setting  the  switch  to  the  desired  baud  rate 

(2)  Repeat  the  break  key  until  an  asterisk  (*)  comes  up;  then 
use  any  keyboard  key  once. 

When  the  baud  rate  is  changed  successfully,  the  message  “SPEED 
CHANGE  COMPLETE"  should  appear. 


The  software  handles  outputs  for  several  different  kinds  of  termi- 


nals. Keying  in  one  of  the  following  parameters  after  typing  CTRL-B  tells 
what  type  of  terminal  is  used: 


parameter 


type  of  terminal 


a 

b 

c 

d 

e 

f 

g 

h 

i 


MINI  (not  avai lable) 

TT33 

T I SILENT  700  (DEFAULT) 
DMEL  DATA  MEDIA 

ADM  LEAR  SIEGLER 

HP45  HEWLETT-PACKARD 
TEC  TEC  INC. 

TEK  TEKTRONIX 

ANN  ANN  ARBOR 


For  example,  by  typing  (ctrl-b)e,  the  message  @ADM  will  appear  on  the  screen; 
then,  whenever  the  output  of  a program  exceeds  2A  lines,  it  will  stop,  ring 
the  bell,  and  wait  until  receiving  ctrl-y  to  continue. 

To  users  planning  on  duplication  of  our  RJET,  please 
take  note  of  the  following  concerning  the  hardware  and  software  of  this  network: 


(1)  The  PD P 1I/0A  is  physically  noisy  with  two  whining  blowers; 
therefore,  the  location  must  be  taken  into  consideration. 

(2)  The  2400  baud  rate  is  recommended  when  running  TEKPLOT  on  a 
Tektronix  terminal.  If  A800  baud  rate  is  used,  set  CPS  to  960  before  doing 
any  plotting. 


8 channel s 


t 

to  LBL 


POP  11/50 


POP  1 1/3** 


DMC I 1 -AR 


DMC11-DA 


BeU  208A 
DATAPHONE 


CDC  7600 
(A  machine) 


CRAY 

(C  machine) 


A800  Baud 
"telephone 
1 ine 


FIG.  1 


Bell  208A 
DATAPHONE 


DMC1 1-DA 
DMC I 1 -AR 

PDP  11/04  with 
1 6k  words  MOS 

DZ-11  DZ-11 


8 channels 


5 Lear  Siegler 
ADM  terminals 


Texas  Instruments 
Silent  700 


Tektronix  4013 


PDP 1 1 /DMA 
Control ler 
Mode  I 121 


Versatec 
Matrix 
D1 100A 
(printer) 


Section  IV 

PLASMA  SIMULATION  TEXT 


Chapters  1 - 7 of  Part  I , Primer  and  Chapters  8 - 10  of  Part  I I , 
Theory  have  been  revised.  Copies  of  these  will  be  made  to  fill  requests 
and  for  use  in  courses  (e.g.  here  and  at  Stanford). 


Section  V 

SUMMARY  OF  REPORTS,  TALKS,  PUBLICATIONS  IN  PAST  QUARTER 


None  thi s quarter. 


VISITOR:  Dr.  Liu  Chen  from  PPPL  visited  with  us  during  May,  completing: 

(a)  "Theory  of  Short  Wavelength  Drift  Tearing  Instabilities" 

(with  P.  H.  Rutherford) 

(b)  "Effects  of  Curvature  on  Drift  Tearing  Modes"  (with  H.  L.  Berk). 


ERRATA 


The  code  ROOTS  was  noted  to  work  incorrectly  for  ring  distribu- 
tions with  the  local  approximation.  M.  J.  Gerver  has  found  an  error  which, 
when  removed,  Y.  Matsuda  found  gave  the  same  results  as  the  code  by  N.  Lind- 


gren.  Change 


001879  SUM2  =>  SUM2*  ( 1 . + OMEGA*EPSl/KYAI ) 


001879  SUM2  =■  SUM2*(KYAI  + 0MEGA*EPSl) 

'm.  H.  Gerver  et  al . , "Normal  Modes  of  a Loss  Cone  Plasma  Slab  with  Steep  Density 
Gradient",  Memo.  No.  ERL-M602,  2k  Sept.  1976,  Appendix  I. 
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